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1 Introduction 

The magnetic field permeating the solar at- 
mosphere governs much of the structure, mor- 
phology, brightness and dynamics observed 
on the Sun. The magnetic field, especially 
in active regions, is thought to provide the 
power for energetic events in the solar corona, 
such as solar flares and Coronal Mass Ejec- 
tions (CME) and is believed to energize the 
hot coronal plasma seen in EUV or X-rays. 
The question remains what specific aspect of 
the magnetic flux governs the observed vari- 
ability. To directly understand the role of the 
magnetic field in energizing the solar corona, 
it is necessary to measure the free magnetic 
energy available in active regions. The grant 
now expiring has demonstrated a new and 
valuable technique for observing the magnetic 
free energy in active regions as a function of 
time. 

The magnetic energy of an active region 
is derived from the magnetic virial theorem 
(e.g. Molodenskii, 1969; Aly, 1984; Low, 
1984; Aly, 1989) which gives the total (force- 
free) magnetic energy as 

E ff = f {xB x + yB y )B z dxdy , (1) 

where B x , B y , and B z are the horizontal and 
vertical components of the magnetic field, x 
and y are the horizontal coordinates on the 
Sun (e.g. Gary and Hagyard, 1990), and zo is 
the lower boundary. When assumptions that 
(a) the magnetic field is force-free, (b) that 
all the magnetic field crossing the chromo- 
spheric boundary is measured, and (c) that 
the coronal magnetic field falls to zero some- 
where outside the active region are imposed, 
the virial theorem yields the total magnetic 
energy in the coronal volume threaded by 
the field lines crossing the lower boundary. 
If the lower boundary is measured at an at- 


mospheric level where the magnetic field is 
force-free, then the magnetic virial theorem 
may be invoked to measure the total magnetic 
energy above active regions. Subtracting the 
energy of the equivalent potential field (with 
the same vertical field as that observed) from 
the total magnetic energy yields the free en- 
ergy in the coronal volume, i.e., the energy 
that is available to power solar activity. 

At first glance it appears that the mag- 
netic energy will depend on the coordinate 
system chosen, since equation (1) contains 
the coordinates x and y. However, if the mag- 
netic field is force-free on the boundary, the 
horizontal Lorentz forces, given by 

F x = [ B X B Z dxdy 

47T J z—zq 

F y = f B y B z dxdy , (2) 

y 47T Jz=Z0 

are zero. Hence, equation (1) is independent 
of the coordinate system and valid when the 
field is force-free. 

Measurements of active-region photo- 
spheric vector magnetic fields have become 
routine. Quantities derived from the pho- 
tospheric magnetic field have been used for 
years to characterize solar active regions and 
their propensity for solar activity. However, 
the measurement of the active-region photo- 
spheric vector magnetic field limits our phys- 
ical understanding of the solar atmosphere in 
a fundamental way: the photosphere is still 
forced, i.e., the plasma /3 = P/8ttB > 1 and 
J x B can be non-negligible. Hence, the use 
of the photospheric measurement to under- 
stand the coronal magnetic structure is com- 
promised to an uncertain degree. 

A method for directly measuring the 
magnetic free energy in active regions was pi- 
oneered by Metcalf et a 1. (1995), based upon 
the chromospheric vector magnetic field de- 
rived from observations of the Nal D-line at 
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Figure 1: The magnetic free energy is measurable only with chromospheric data: only the 
data point highest in the atmosphere (closest to line center wavelength) samples the field in 
a force-free region. Here the virial theorem gives a valid measure of the magnetic energv in 
AR 7216. 


5896 A. This magnetically sensitive spectral 
line is formed high enough so that the mag- 
netic field is force-free (J x B « 0) but low 
enough so that the Zeeman polarization sig- 
nals are sufficiently strong for a reliable mea- 
surement. In the initial Metcalf et a 1. (1995) 
study, the Stokes polarization spectra across 
the Nal line were recorded with the Haleakala 
Stokes Polarimeter (HSP) at Mees Solar Ob- 
servatory, Hawai‘i as the instrument rastered 
across NOAA Active Region (AR) 7216. 

Metcalf et a 1. (1995) determined the 
height at which the measured field becomes 
force-free. When the magnetic field is not 
force-free, Lorentz forces (Eq. 2) are non- 
negligible and the virial theorem breaks 
down. Figure 1 demonstrates this effect: 


magnetograms constructed from the polariza- 
tion signal at wavelengths farther from line 
center ( i.e ., deeper in the solar atmosphere) 
result in error bars that become exceedingly 
large and a computed total energy which is 
less than the energy of the equivalent poten- 
tial field. This unphysical state is directly due 
to the breakdown of the assumptions listed 
above, i.e., an indication that the atmosphere 
in these deeper layers is not force-free. 

Hence, we must observe the magnetic 
field in the chromosphere, high enough above 
the photosphere so that the field is force- 
free, if we are to routinely apply the virial 
theorem. Figure 1 demonstrates that, below 
AA ~ 80 mA from Nal line-center, the poten- 
tial field energy drops below the total energy, 
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Figure 2: The IVM chromospheric line-of-sight magnetic field in AR 8299 overlaid on a 
continuum image. Red contours show fields directed towards the observer and blue contours 
show fields directed away from the observer (100, 200, 400, 800 G). Yellow contours show the 
latitude and CMD; the axes are labeled in arc seconds. The dark triangles and the apparently 
strong magnetic field which appear at the center of the extreme upper and lower edges are 
from fiducial marks in the IVM and are neither real nor included in our calculations. 


and by AA f»68 mA the chromospheric mag- 
netic field fulfills the force-free criterion and 
provides a reliable measure of the magnetic 
energy using the virial theorem. 

2 Recent Advances 

2.1 Multi-wavelength observa- 
tions 

Under this grant, the IVM and Mees Obser- 
vatory conducted three separate campaigns 


to observe in the chromospheric Na-D line. 
Table 1 lists the dates of these campaigns 
and Figure 2 shows a single example from this 
dataset. 

To advance the science goals of this 
grant, we undertook these campaigns to ob- 
serve chromospheric magnetograms with the 
IVM in it’s old configuration. From 7-Aug- 
1998 through 4-Sep-1998, a total of 29 ac- 
tive regions were surveyed and 17 movies se- 
quences recorded. From 19-Dec-2000 through 
l-Feb-2001, a total of 41 active regions were 
surveyed, and 12 movie sequences recorded. 
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Data from these campaigns is under active 
analysis. 

The HI H-alpha line was observed dur- 
ing 26-Jul-2001 through 09-Aug-2001. This 
was a first attempt to make IVM magne- 
tograms in the Ha line. Despite a number 
of setup problems, quick-look analysis shows 
that valid signal was detected in all Stokes 
parameters in many active regions. Analy- 
sis of these data will be conducted to assess 
the utility of H-alpha observations for char- 
acterizing the chromospheric field at a higher 
altitude than is seen in Nal -D. 

2.2 IVM Hardware Modifica- 
tions 

2.2.1 Tunable filter development 

The IVM air-gap Fabry-Perot etalon is a high 
quality spectrometer. With a finesse of 50, 
it is one of the best operational etalons in 
regular use for solar observations. Because 
it is designed to provide high spectral resolu- 
tion FWHM of 70 mA, the free spectral range 
(the separation between etalon orders) is 3.5 
Angstroms. Therefore, we need a blocking fil- 
ter that has a FWHM of 3 Angstroms or less. 
Such narrow filters are available, but limit the 
range of spectral scans within a single line, 
reduce overall system transmission, and are 
sensitive to alignment and temperature vari- 
ation. It was desirable therefore to develop 
a tunable blocker that could be used in con- 
junction with the etalon, and would permit 
rapid switching among widely separated spec- 
tral lines. 

As originally designed, the IVM was to 
use an acousto-optic tunable filter. That 
technology did not mature, and we had to 
fall back to fixed interference filters as block- 
ers. The time to open the system, remove 


the blocker, install a new one, close and sta- 
bilize the system is one day. Therefore, com- 
parisons between photospheric and chromo- 
spheric magnetograms and/or structure im- 
ages are not possible. 

Under this grant, we pursued a new tech- 
nology: liquid crystal tunable etalons. This 
technology was initially developed for use 
in aeronomy, both ground and space-based 
(Schneller et al., 1996) The potential advan- 
tages of these etalons are many: 

• Wide angle of acceptance. The index of 
refraction of the liquid crystal is in the 
range 1.5 to 1.8, increasing the angle 
compared to a gap etalon. 

• Rapidly tunable at low voltage. 

• Stable alignment. Gap etalons regu- 
larly lose parallelism and must be re- 
aligned. Liquid crystal gaps are set by 
mechanical spacers. 

• Low cost. 

• Compact design. Gap etalons require 
significant structure around the clear 
aperture to hold the drivers and sensors 
for the tuning system. 

The liquid crystal etalons are produced 
by Scientific Solutions, Inc. Initial discus- 
sions with the vendor demonstrated that 
many of the basic problems unique to our ap- 
plication were tractable. In aeronomy appli- 
cations, overall transmission is critical, there- 
fore reflector coatings are held to modest re- 
flectance in order to minimize absorbance. 
Typical etalons then have a finesse of 15, 
reflector limited. The surface flatness and 
parallelism of the substrates appeared not to 
limit the finesse. 

For the IVM, we decided to attempt to 
accomplish both rapid switching among lines 
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Dates 


Spectral line Comments 


03-Jul-1996 - 30-Jul-1996 Nal 

07-Aug-1998 - 04-Sep-1998 Nal 

19-Dec-2000 - 01-Feb-2001 Nal 

26-Jul-2001 - 09- Aug-2001 Ha 


29 Active Regions 17 Movie sequences 
41 Active Regions 12 Movie sequences 


Table 1: Nal Data were obtained in three separate observing campaigns. 


The Magnetic Free Energy in Active Regions, Thomas R. Metcalf and Donald L. Mickey 
1999), American Astronomical Society Meeting 194, 9407 

A Comparison of the Active Region Magnetic Field in the Photosphere and Chromosphere., 
Thomas R. Metcalf and K. D. Leka 2001, American Geophysical Union, Spring Meeting 
2001, abstract #SP41B-07 441B07 

A Comparison of the Active Region Magnetic Field in the Photosphere and Chromosphere, 
Thomas R. Metcalf and K. D. Leka 2001, American Astronomical Society, Solar Physics 
Division Meeting 

The Magnetic Free Energy and a CME in Active Region 8299, Thomas R. Metcalf, Donald 
L. Mickey, Barry J. LaBonte, and Leigh Ann Ryder, Astrophys. J., submitted, 2001. 

Active Region Magnetic Structure Observed in the Photosphere and Chromosphere, K. D. 
Leka k Thomas R. Metcalf 2001, Astrophys. J., in preparation 


Table 2: Publications acknowledging this grant. 


and broader tuning within a spectral line. 
This requirement can be met with a pair of 
etalons, each having finesse 30. Models of 
the solar spectrum times the etalon transmis- 
sion were computed by us and by the vendor 
as a check on our mutual understanding of 
the materials properties and the etalon per- 
formance. After many iterations, a design for 
a pair of etalons used in vernier mode was 
settled upon as optimal. This would permit 
scanning to any wavelength between 580 and 
660 nm, with a parasitic light level of a few %, 
no more than the present IVM. By using the 
availability of standard spacers and offsetting 
the index of refraction of the liquid crystal, a 
wide choice of vernier ratios could be tried. 
The optimal design etalon pair had spacer 
thicknesses of 10 and 15 microns, and indices 
of 1.79 and 1.58 respectively, for a vernier ra- 


tio of 1:1.32. 

The vendor worked out a design for a 
higher reflectance coating that would permit 
higher etalon finesse. The coating was de- 
posited on a substrate, its reflectance verified, 
and the distortion of the substrate measured. 
The substrates for the etalons were then fig- 
ured with a compensating distortion. When 
they were coated, they then are flat. Both 
the reflectance and surface figure appeared to 
meet the requirements. 

The etalons were assembled and tested. 
One had a measured finesse of 18, and the 
other a finesse of 10. The overall transmit- 
tance was lower than expected. The cause of 
the finesse failure is not certain, but may be 
related to the two additional steps that oc- 
curred after the reflectors were tested: depo- 
sition of the transparent electrodes, and wip- 
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ing of the surface to align the liquid crystal 
molecules. The cause of the transmission fail- 
ure is likely higher than expected absorbance 
in the coatings, an effect difficult to charac- 
terize by direct measurement. 

After discussion of alternatives, it was 
decided that this technology is not yet ca- 
pable of the higher requirements that the 
IVM application places upon it. It remains 
a technology of high potential, and new ma- 
terials handling techniques, such as atomic- 
beam alignment of the LC molecules (Chaud- 
hari et al., 2001) in place of mechanical rub- 
bing may improve the performance for appli- 
cations similar to the IVM. 

2.2.2 Filter Wheel 

The fail-back plan against a failure of the LC 
etalons was to install a filter wheel in the 
IVM Data Camera beam. The original IVM 
blocker filters were fabricated in the early 
1990’s; nearly a decade of progress in that 
technology made it possible to procure filters 
with higher transmission and steeper cutoff of 
offband light. However, limited space is avail- 
able in the IVM instrument where the wheel 
must be located. To accommodate the wheel, 
a new sub-plate had to be fabricated. 

The mechanical parts for the wheel are 
complete. The motor, drive, and filter ther- 
mal control components are procured. The 
IVM observing software was designed with 
the components for a tunable filter, which 
will be used to sequence the wheel with the 
observing program. We expect to have the 
wheel test run within a month, as the last 
elements are integrated. At a later time, we 
will change the final camera lens and add a 
lens adjacent to the wheel to render the beam 
telecentric, improving the uniformity of the 
blocker filters over the FOV. 


2.3 Enhancement to the IVM 
Data Reduction Software 

At the start of this grant, the IVM data re- 
duction software was not capable of analyz- 
ing the Nal data. The first analysis task that 
faced us was to upgrade the IVM software to 
allow the analysis of the data listed in Ta- 
ble 1. The software was originally designed 
to work with the relatively narrow Fel lines 
and the broad Nal line required substantial 
modifications to the analysis software. The 
enhancement of the software is complete and 
the data reduction software is now capable of 
routine analysis of Nal data. 

2.4 Observations of the 3-D 
Structure of Sunspots 

Advances in the inversion procedures of 
Stokes spectra have enabled photospheric- 
line spectropolarimetry to reveal the depth- 
dependence of a number of quantities of in- 
terest. Westendorp Plaza et al. (2001a, 
2001b) have demonstrated techniques which, 
when applied to the 630.15,.25nm Fe-I pho- 
tospheric lines, result in depth-variation in- 
formation of the vector magnetic field, tem- 
perature, and velocity stratification between 
—2.5 < logr 5 < 0. This is roughly a vertical 
span of 300km from deep in the photosphere 
to just below the chromosphere. The results 
of this approach are promising but still lim- 
ited to a narrow range of atmospheric heights. 
Additionally, these authors focus squarely on 
the variations observed in the penumbra of a 
small sunspot. 

A first direct comparison of vector mag- 
netic fields in the photosphere and chro- 
mosphere with near-cotemporal observations 
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arcseconds 


Figure 3: Continuum image of NOAA AR8299 from the Imaging Vector Magnetograph ob- 
tained near the Na-I 589.6nm line. Contours show the sunspots observed with the Advanced 
Stokes Polarimeter in the continuum near the Fe-I 630nm lines. 


used ASP and IVM data from NOAA 
AR8299 (N19 W62) on 18 August 1998 
(Fig. 3). The Advanced Stokes Polarime- 
ter scanned this active region from 13:44— 
13:56UT in the photospheric Fe-I lines de- 
scribed above. Observations with the Imag- 
ing Vector Magnetograph scanning the chro- 
mospheric Na-D line began at 17:23 UT and 
continued for many hours with a 3-5 min. 
cadence. For the initial analysis, we in- 
creased the IVM signal-to-noise ratio by aver- 
aging the spatially-coaligned spectra of eleven 
datasets from 17:23-17:53, all of which had 
good seeing. 

The preliminary comparison of these two 
datasets (to be submitted for publication in 
Autumn 2001) was work funded primarily by 
NASA through NASG-00031 and this grant. 


We present the main results here. First and 
foremost, a comparison of the magnetic flux 
observed at the two layers is shown in Figs. 
4, 5 and 7. 

The observed changes between the pho- 
tospheric and chromospheric vector magne- 
tograms include a weakening of the total mag- 
netic flux at the higher layers (Figs. 5, 7). 
This is hardly surprising, but the magnitude 
of the weakening is. Note that the compar- 
isons here are limited to quantities derived 
without performing the ambiguity-resolution 
of the transverse component of the magnetic 
flux. The reason for this limitation is that, for 
observations close to the limb, the ambiguity 
resolution ( e.g Metcalf, 1994) requires mix- 
ing between the transverse and line-of-sight 
fields to derive the horizontal and vertical 
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Photospheric Magnetogram Chromospheric Magnetogram 



Figure 4: Magnetic Flux maps showing the line-of-sight (contours: = positive, at 

50,100,200,400,800, 1600G levels) and transverse (line-segments), magnetic flux components. 
Top: photosphere, binned to 1" pixels, maximum B t =2500G. Bottom: chromosphere, maxi- 
mum B t =1700G. 


fields and the ambiguity resolution becomes 
less reliable. This is obviously a limitation 
of this dataset (and a strong motivation for 
additional coordinated observations). 

Additionally, the direct comparison of 
the two datasets clearly shows the super- 
penumbral structure of the sunspots. Wes- 
tendorp Plaza et a 1. (2001a) detected the 
sunspot canopy as an increase in the magnetic 
field strength with height (through the upper 
photosphere) beyond the visible penumbral 
boundary. We detect the magnetic canopy or 
superpenumbra directly in the chromosphere, 
as a positive fractional change in total flux 
(Fig. 7) occurring outside the photospheric 
penumbral boundary (Fig. 8). 

2.4.1 Heights of Formation and Heights 
of Detection 

At disk center, the Fe-I lines are formed over 
the range 250-350km above r 5 oonm = 1 (Bruls 


et a 1., 1991), and, in the VAL-F model atmo- 
sphere, the Na-D line is formed at approx- 
imately 2 = 1000km above r 5 oonm = 1 at 
68mA from line-center in the quiet sun (e.g. 
Metcalf et al., 1995). Hence, from quiet Sun 
models at disk center, a first estimation of 
A 2 , the difference in the height of formation 
of Fe-I and Na-D, is of order 700 km. A cor- 
rection to this must be made for the obser- 
vations presented here, in that the observing 
angle fj, = 0.47. Spectral features which form 
at optical depth unity form at r vert = 1/0.47 
when the atmosphere is viewed at an angle 
corresponding to // = 0.47. This is typically 
less than one density scale height higher than 
^vert ~ 1 (assuming an exponentially strati- 
fied atmosphere). In a non-magnetic atmo- 
sphere this would be only 100 km. However 
in a sunspot’s magnetic atmosphere the vari- 
ation in A z could be significantly greater. 

To better grasp the A z question we per- 
formed a simple test. The active region was 
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Figure 5: Comparison of the photospheric horizontal axis and chromospheric (vertical axis) 
line-of-sight (left) and transverse (right) magnetic flux. Error bars are omitted for clarity, 
however points below the noise levels are not plotted. 


very close to a potential configuration, with 
a region-averaged twist parameter a rm AR ~ 
—0.008 Mm - 1. Extrapolations in height 
from the photospheric flux were performed, 
and compared to the observed chromospheric 
flux. The goal was to find the height at which 
they best agreed (Fig. 9). The results are 
surprising: the best agreement comes with a 
A z =2-2. 5Mm. 

2.4.2 Magnetic scale heights 

The logical extension from these initial re- 
sults is to investigate dB/dz and the im- 
plied magnetic scale heights. Using the 
mean Az derived from the quiet-sun mod- 
els (0.7 Mm) and the extrapolations (2-2.5 
Mm), and allowing for the possible range 
implied by their substantial difference, we 
show in Fig. 10 the inferred dB/dz across 
different sunspot components. We find val- 
ues which are generally lower than the « 2 
G/km derived from deep-photosphere/upper- 
photosphere observations (Westendorp Plaza 
et ai.2001a; Pahlke & Wiehr 1990) but 


which agree with previous observations of 
the photosphere/chromosphere values, gener- 
ally « 0.5 ± 0.2 G/km with lower values in 
the penumbra! than umbral regions (Ruedi, 
Solanki, & Livingston 1995a; Ruedi et a i. 
1995b). 

Inverting this quantity (dB/dz) and nor- 
malizing by the photospheric magnetic flux 
gives us the implied magnetic scale heights 
(Fig. 11). The most obvious feature in the 
scale height is the uniformity in the umbra, 
demonstrating that the umbra is a monolithic 
structure filled with magnetic field. Such uni- 
formity is expected given the large magnetic 
fill-factor in the photosphere and the consis- 
tency of the field strengths across sunspot 
umbrae. The second feature in the implied 
scale heights is a dramatic change once the 
umbral/penumbral boundary is crossed, from 
almost a single value in the umbra to a wide 
range covering almost an order of magnitude 
in the penumbra and plage. 

The puzzling result is that no mat- 
ter what end of the range of possible 
Az used, the scale height is far differ- 
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Figure 6: A comparison of photospheric and chromospheric Stokes profiles observed with the 
ASP and the IVM. The broad chromospheric Nal line is particularly well suited to a weak 
field analysis. 


ent from what is expected. Generally 
the magnetic scale height, especially in the 
umbra, would be of the same order as the 
size-scale of the sunspot umbrae itself; in this 
case, at least lOMm or 10 4 km. With a sim- 
ple dipole model, the sunspot scale height 
increases linearly with height from a value 
roughly half the size of the umbra at the pho- 
tosphere to a value consistent with the size of 
the umbra at a height which is also the same 
as the umbral diameter. Hence, there is a dis- 
crepancy of at least a factor of two between 
our results and the size of the sunspot itself 
(the sunspot magnetic field is decreasing with 
height more rapidly than expected), unless we 
neglect the atmospheric models and assume 
that the height we derive from the potential 
extrapolations is correct. 


2.4.3 The Puzzle 

The conundrum is simple: sunspot mod- 

els which do not explicitly include magnetic 
fields indicate that chromospheric Na-D emis- 
sion should come from a certain atmospheric 
height, while extrapolations from the photo- 
sphere lead to a much larger height. 

The modus operandi until now has 
been to employ photospheric magnetic flux 
data as a boundary condition for coronal- 
field extrapolations. This approach explic- 
itly ignores the transition from a plasma- 
dominated to magnetic-dominated regime 
(0 = Pg»a/ Pmag > 1 - ► fd < 1). Pre- 
liminary evidence from our dataset points 
toward less of a discrepancy between coro- 
nal field extrapolations and coronal X-ray 
and EUV observations when chromospheric 
(force-free), rather than photospheric (non- 
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Figure 7: Percent change in total magnetic flux B — \j Bf + B\ plotted as a function of 
continuum intensity (normalized by the quiet-sun mean intensity). The error bars reflect 
the uncertainties in both the photospheric and chromospheric flux measurements, and only 
points above the noise are considered. The (negative) fractional change in the umbra is 
a slight function of continuum intensity but consistent with a constant value. This is in 
contrast with the wide range of fractional change in the penumbra and plage. The points 
with positive fractional change, i.e., with chromospheric flux greater than photospheric, occur 
only in the edge of the penumbra and in the plage regions, according to continuum intensity. 
As shown in Fig. 8, these points represent the direct detection of the superpenumbra, or 
canopy-like extension of magnetic field above and beyond the photospheric penumbra. 


force-free), magnetic field measurements are 
used as boundary conditions for interpreting 
coronal (probably still force-free) conditions 
(Figure 12). Hence, the “problem child” is 
the photosphere/chromosphere transition. 

Further, the interpretation of data from 
different atmospheric heights has routinely 
assumed single- values or single-ranges for the 
formation heights of a spectral line across 
very disparate solar structures. Observations 
are now forcing us to reconsider the simplistic 
assumptions and approaches of the past. 

Sunspots are one building block of the 
solar magnetic cycle, exerting their influence 


over an enourmous atmospheric range. These 
bundles of strong magnetic fields influence so- 
lar accoustic signals below the photosphere 
and push aside extant magnetic structures 
into the corona. Active regions carry free 
magnetic energy which may power solar ac- 
tivity such as flares and Coronal Mass Ejec- 
tions. The work funded by this grant has now 
forced us to confront some glaring inconsis- 
tencies in our understanding of basic sunspot 
structure. 

2.5 Observations of the Mag- 
netic Free Energy 
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Delta B (Chromosphere — Photosphere) 



Figure 8: The super-penumbra as seen in a difference image between chromospheric and 
photospheric magnetic flux. Dark areas have negative fractional flux change (see Fig. 7), 
light areas have positive flux change and are highlighted by dark contour. Lighter contours 
indicate photospheric umbral and penumbral boundaries. Both the structure of the fractional 
flux variation and the location of the superpenumbra are visible in this example. 


2.5.1 Confirmation of the Technique 

Metcalf et al. (1995) measured the magnetic 
free energy with a single observation of a sin- 
gle active region, demonstrating the feasibil- 
ity of this type of observation. While the 
technique successfully measured the free en- 
ergy in AR 7216 at 10 33 ergs (Figure 1), the 
vector magnetogram required nearly 5 hours 
of observations, making the HSP an imprac- 
tical instrument for studying variations of the 
magnetic free energy. 

With funds from this NASA grant, the 
U. Hawai‘i/Mees Solar Imaging Vector Mag- 
netograph (“IVM”) was used to observe 
the chromospheric vector magnetic field in 
AR 8299 (Figure 2). The IVM (Mickey 
et al. 1996; LaBonte et al. 1999) has 


a symmetric design with near-normal reflec- 
tions and a helium-filled telescope to mini- 
mize instrumental polarization and internal 
seeing. A four-frame polarization-modulation 
sequence is used; these mixed polarization 
states are sampled by means of a Fabry Perot 
etalon at 30 wavelength positions with 70mA 
FWHM spectral resolution across the mag- 
netically sensitive line of choice; the default 
is the photospheric Fel 6302. 5A line com- 
monly used for photospheric vector magne- 
tograms. By changing the prefilter, the IVM 
can obtain data in other spectral lines, includ- 
ing the Na I chromospheric line. The IVM 
was designed for active-region observations; 
it records spectropolarization images with 
512 x 512 0.55" pixels covering over 4' x 4'. A 
single spectrum-scan is obtained in under two 
minutes; to increase the signal/noise ratio, it 
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Figure 9: Comparison of the chromospheric (horizontal axis) B Jong (left) and B.trans (right) 
magnetic flux to extrapolated flux (using the photospheric flux as the boundary). Bottom 
panels, the correlation coefficients for line-of-sight (left) and transverse (right) fluxes between 
the two, as a function of A z between the two. Best fits, weighting the lower-noise higher- 
correlation B i more than B t , is Az =2-2. 5Mm. This contrasts with the expected « 0.7Mm 
from non-magnetic stratified atmospheres (see text). 


is common to average the spectra from two 
or more coaligned data sets before inversion 
to the magnetic flux vector. 

On 1998 August 11, the IVM observed 
AR 8299 in the Nal spectral line for the 
Whole Sun Month campaign (Figures 2,13). 
From these chromospheric magnetograms, 
the total magnetic energy was computed as a 
function of time. The free energy is the total 
energy less the energy of the equivalent po- 
tential (current-free) magnetic field; the po- 
tential field energy, computed using the ob- 
served vertical field as a boundary condition 
(Fig. 14), was nearly constant at 10 32 ergs. 

The initial results confirmed the tech- 
nique: first and foremost, the magnitude of 
the free energy in AR 8299 is consistent with 
the work of Metcalf et al. (1995) which used 
a very different instrument. The primary dif- 
ference between this observation and the ear- 


lier work is the considerably smaller value 
of the potential magnetic energy found in 
AR 8299: the total energy is nearly an or- 
der of magnitude larger than the potential en- 
ergy. The magnetic energy is about a factor 
of two lower than the open field limit given by 
Aly (1984), indicating that the large observed 
free energy is physically possible (dashed line 
in Figure 14). 

The secondary difference is that the er- 
ror bars for AR 8299 are somewhat larger 
than for AR 7216; this is due to some of the 
magnetic flux leaving the field-of-view and 
connecting to nearby AR 8297 (Figure 13). 
Recalculating the total energy for many dif- 
ferent origins of the coordinate system will 
yield the same result unless the field is not 
force-free and/or significant flux is leaving 
the field-of-view. Any statistically significant 
variation of the energy with the coordinate 
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dB/dz (and RANGE implied by dz uncertainty) 



Figure 10: Inferred dB/dz from the photospheric and chromospheric total magnetic flux, 
as a function of Ic/Iqs. Points show the mean, “error bars” indicate the possible range of 
dB/dz given the possible range of formation heights for both the Fe and Na lines. 


system will quantify these systematic effects. 
For AR 8299, the magnetograms are derived 
from AA=68mA from line-center where the 
field is believed to be force-free. Additionally, 
a careful analysis of the flux loss to AR8297 
using KPNO data indicates that the flux loss 
is acceptable and the virial theorem applies. 

2.5.2 Variation in the Magnetic Free 
Energy 

As mentioned above, the IVM has the capa- 
bility of obtaining a full vector magnetogram 
every few minutes, enabling not only the mea- 
surement of the magnetic free energy but also 
its variation in time: Figure 14 shows the to- 
tal magnetic energy as a function of time for 
AR 8299. The free energy peaks at about 
1.5 ±0.5 x 10 33 ergs, but drops to a level con- 
sistent with zero for almost an hour late in the 
observation before gradually increasing to the 
previous level. As the observation ended, the 


total magnetic energy approached the open 
field limit, suggesting that the field was in a 
substantially more open configuration com- 
pared to the field before the dip. 

This intriguing observation will be dis- 
cussed in some detail below since it is a good 
example of the type of observation that we 
have obtained with funds from this grant. 
The obvious question is, what led to this dra- 
matic change in the magnetic free energy? 

2.5.3 An Invisible CME? 

The remarkable drop in the magnetic free en- 
ergy observed between 20:30 and 21:30 UT 
was accompanied by almost insignificant soft 
X-ray activity on the Sun (Figure 15). There 
were no radio bursts reported during this 
time interval, the Mauna Loa coronagraph 
was not observing, and CGRO/BATSE saw 
no hard X-ray bursts. Unfortunately, no 
SOHO data are available and TRACE was 





Figure 11: Similar to Fig. 10 but showing the implied magnetic scale height. The constancy 
of this scale height is remarkable over the umbra; the abrupt change in scale heights between 
umbra/penumbra is striking, although we use a constant range of Az for all points regardless 
of their place in the sunspot. 


not observing this active region. Hence, we 
have little information on what may have oc- 
curred in the corona of this active region dur- 
ing this time. 

However, Figure 16 shows SXT compos- 
ite images before and after the dip in the 
magnetic free energy (no SXT data are avail- 
able for AR 8299 between 19:38 UT and 22:34 
UT due to calibrations). There is a signif- 
icant change in the X-ray structure of the 
active region at some time between the im- 
ages in Figure 16 which includes the energy 
dip. The most obvious change is the appear- 
ance of a bright cusp connecting AR 8299 and 
AR 8297, seen in the center of the SXT ROI 


difference image. Unfortunately, we do not 
know exactly when this connection bright- 
ened. We also note that there was a faint soft 
X-ray sigmoid structure connecting AR 8299 
with AR 8297 prior to the dip in the magnetic 
energy, though this structure does not appear 
to change between 19:38 and 22:34 UT. 

The existence of an X-ray sigmoid, the 
development of a cusp structure, as well as 
the lack of a GOES event, are all signatures 
often associated with gradual, halo CMEs 
(e.g. Webb et a 1. 1998; Hudson et a 1. 1998; 
Canfield, Hudson &: Pevtsov, 2000; Sheeley 
et a 1. 1999). Hence, the SXT data pro- 

vide strong support for the occurrence of a 
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Figure 12: The magnetic field extrapolation from the chromosphere (top) matches the SXT 
data better than the extrapolations from the photosphere (bottom). 


gradual CME. A CME without any remark- 
able effects in the low corona is not unheard 
of: Webb et a 1. (1998) report such an event 
from 1997 January 6, in which a geoeffective 
CME was associated with remarkably weak 
coronal activity (a GOES Al event) and a 
combination of brightening and dimming in 
SXT images. This scenario is very reminis- 
cent of the 11 August 1998 AR 8299 observa- 
tion presented here. Without any data from 
SOHO/LASCO, we cannot be certain that 


such a CME occurred in AR 8299, however 
the circumstantial evidence suggests that a 
CME was launched. 

The observation presented above is, to 
our knowledge, the first direct measurement 
of the time development of the magnetic free 
energy in an active region. Relating the ob- 
served changes in the free energy with activ- 
ity on the Sun, we may have observed the 
direct energization of a gradual CME for the 
first time. The dramatic drop in free mag- 
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Figure 13: The line-of-sight magnetic field overlayed on an SXT image (negative image, 
logarithmically scaled). The gray square shows the IVM field-of-view. The magnetic field is 
the same as Figure 2 and the axes are labeled in arc seconds. AR8299 is in the center of the 
image and AR8297 is the upper right region. 


netic energy observed in AR 8299 amounts to 
approximately 10 33 ergs. This is in an enor- 
mous amount of energy given that there was 
no indication of flaring or other activity at 
the time of the energy dip. Where did the 
energy go? 

While it is unlikely that a gradual CME 
could account for the full 10 33 ergs (e.g. Anti- 
ochos et a i. 1999), the topological change as- 
sociated with a CME might redistribute the 
free energy out of the active region we ob- 
served. This hypothesis is supported by the 
rapid restoration of the free energy starting at 
21:00 UT. The increase in free energy from 
21:00 to 22:00 amounts to about 10 33 ergs, 


implying a power input of 3 x 10 29 ergs s 1 . 
While possible, it is unlikely that this energy 
build-up could come from photospheric mo- 
tions over such a short timescale, since pho- 
tospheric motions are typically only 1 km s -1 
(e.g. equation (27) in Aly, 1989). Further, 
there is no observational evidence for emerg- 
ing flux which might stress the coronal field 
more quickly. Hence, the redistribution of 
free energy within the AR 8299/8297 sys- 
tem, combined with a gradual CME, is a 
more likely explanation of the variation in 
the free magnetic energy. New observations 
which can be combined with Yohkoh/SXT 
and SOHO/LASCO data are required to fully 
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Figure 14: The total magnetic energy (ergs) above the chromosphere in AR 8299 (solid line) 
dips to a value consistent with zero around the time of a CME. The vertical error bars 
indicate both statistical and systematic errors. The horizontal error bars indicate the time 
interval over which each observation was averaged. The dotted line shows the energy of the 
equivalent potential field and the dashed line shows the equivalent open field energy. 


understand the remarkable changes possible 
in solar active regions. 

2.5.4 Where’s the Flare? 

The storage, loss, and subsequent gain of over 
10 33 ergs in this active region begs the ques- 
tion, why did it not produce a flare? There is 
sufficient free energy to power even a large so- 
lar flare, but energy storage alone is not a suf- 
ficient condition for the production of a flare. 
This observation verifies that, not only must 
energy be available in the magnetic field, but 
a proper field topology and a trigger must 


also be present. Apparently these conditions 
were not fulfilled in AR 8299. The precise 
conditions required for flares are not well un- 
derstood and we expect further studies of the 
changes in magnetic free energy in flaring ac- 
tive regions to shed significant light on this 
subject. 

While the dip in the free energy is the 
most remarkable feature of this dataset, the 
quasi-steady free energy outside this inter- 
val is also interesting. The magnetic energy 
stored in the active region magnetic fields is 
more than sufficient to explain the enhanced 
heating of the active region corona. In fact, 
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Figure 16: Composite, negative Yohkoh/SXT X-ray images of the full disk (logarithmically 
scaled) and AR 8299/8297 (linearly scaled) using the AlMg analysis filter. Left: before the 
energy dip. Middle: after the energy dip; Right: difference image. Around the time of the 
energy dip, a connection between AR 8299 and 8297 brightened dramatically. 


dramatic change in magnetic free energy was 
accompanied by SXT observations suggestive 
of a gradual CME but the data are circum- 
stantial. If a CME did occur in association 
with the drop in magnetic free energy, this 
would constitute the first direct evidence that 
the CME is powered by free energy stored 
in the magnetic field. New observations us- 
ing the Nal D-line in conjunction with di- 
rected, coordinated space-based observations 
from Yohkoh, TRACE, SOHO, and HESSI 
are clearly called for at this time. 

Our next goal will be to determine the 
physical changes which cause variations in 
the magnetic free energy. That is, we have 
seen the effect, but what was the cause ? Ul- 


timately, we hope to make advance predic- 
tions of solar activity by observing variations 
in the available magnetic energy. For a CME 
with little or no signature in X-rays, such as 
the event postulated here or the 1997 Jan- 
uary 6 event studied by Webb et a 1. (1998), 
a dramatic change in the free energy may be 
the only indication, apart from coronagraph 
data, that a CME occurred. 
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Start Time (1 1 -Aug-98 19:00:00) 

Figure 15: The GOES-10 soft X-ray light curves (W m~ 2 ; upper curve: 1-8A, lower curve: 
0.5-4A). There was no significant flaring activity in AR 8299 at the time of the CME and 
the dip in the magnetic free energy. 


the 10 33 ergs of free energy could power the 
observed X-ray emission (10 26 ergs s -1 ) for 
several months, even if it were not replen- 
ished, and it is much larger than a simple 
estimate of the thermal energy (< 10 29 ergs) 
in the active region corona. It has been sug- 
gested that the energetics of the quasi-steady 
enhancement of the heating in the active re- 
gion corona is dominated by the magnetic 
field (Fisher et ad. 1998). In the case of 
AR 8299 we have found that the free energy is 
more than sufficient, by several orders of mag- 
nitude, to counter the X-ray radiative energy 
losses in the active region and to explain the 
enhanced heating of the active region if the 
magnetic energy can be converted to thermal 
energy. 


3 Conclusion and Future Work 

Since the magnetic field in coronal active re- 
gions very likely dominates the energetics of 
coronal activity, measurement of the mag- 
netic free energy available to the corona is an 
important step forward in our understanding 
of active regions. With the support of this 
grant, we have now demonstrated the feasi- 
bility of making temporally resolved measure- 
ments of the magnetic free energy above ac- 
tive regions. We have shown that the chromo- 
spheric magnetic flux vector, measured high 
enough above the photosphere so that the 
magnetic field is force-free, renders the mag- 
netic virial theorem applicable. 

We are now in a position to fully ex- 
ploit this new tool. The observations pre- 
sented herein demonstrate the usefulness of 
the technique, however, the results provide a 
mystery: where did the free energy go? A 
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